Phase diffusion and locking in single-qubit lasers 
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Motivated by recent experiments, which demonstrated lasing and cooling of the electromagnetic 
field in an electrical resonator coupled to a superconducting qubit, we study the phase coherence 
and diffusion of the system in the lasing state. We also discuss phase locking and synchronization 
induced by an additional ac driving of the resonator. We extend earlier work to account for the 
strong qubit-resonator coupling and to include the effects of low-frequency qubit's noise. We show 
that the strong coupling may lead to a double peak structure of the spectrum, while the shape and 
width are determined to the low-frequency noise. 
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In a number of recent experiments (here we only cite 
few examples) superconducting qubits coupled on chip 
to electrical or mechanical resonators displayed quan- 
tum electrodynamic effects and opened the field of "cir- 
cuit QED" [1-9]. By creating a population inversion 
between two charge states in a driven superconduct- 
ing single-electron transistor (SSET) Astafiev et aJ. [8] 
demonstrated lasing behaviour of a microstripline res- 
onator coupled to the qubit. Grajcar et al. [9] coupled 
a driven flux qubit to a low-frequency LC resonator and 
observed both cooling and a trend towards lasing of the 
resonator field. In contrast to usual lasers, where many 
atoms are weakly coupled to the electromagnetic field, in 
single-qubit lasers one artificial atom is coupled strongly 
to the resonator. In addition, solid state qubits are sub- 
ject to decoherence effects. Some of the consequences 
and novel behavior had been analyzed in Refs. [10-15]. 

Even in the lasing state, the coherence of the electro- 
magnetic field is lost due to phase diffusion after a char- 
acteristic time Td [16], an effect which is observable, e.g., 
in the laser spectrum. Phase diffusion can be suppressed 
by injection locking, that is by driving the resonator with 
an additional coherent signal. This fixes the phase differ- 
ence between the laser and driving field to a value which 
depends on the intensity and detuning of the latter. Both 
injection locking and phase diffusion were studied exper- 
imentally for a single-qubit laser in Ref. [8]. As com- 
pared to the spectrum observed in standard (many-atom) 
lasers, the single-qubit laser spectrum is broader and the 
peak is substantially shifted with respect to the natural 
resonator frequency. The maximum photon number in 
the resonator is rather low, raising questions about the 
coherent nature of the amplification in these systems. 

In this work we study the spectral properties of single- 
qubit lasers and explain qualitatively several of the exper- 
imental observations. In Section I we introduce the model 
and describe our approach. Results for static properties 
of single qubit lasers are discussed in Section II. Here 
we focus on the average photon number in the vicinity 



of the lasing transition, which illustrates the differences 
between single- and multi-atom lasers. We show explic- 
itly that single qubit lasers are characterized by a smooth 
transition to the lasing regime and by the absence of a 
sharp lasing threshold. Next, in Section III, we ana- 
lyze the phase diffusion of the resonator field focusing 
on the following issues: (i) we discuss the effects of cor- 
relations between qubit and resonator on the diffusion 
process; (ii) we show how the interplay between strong 
coupling and spontaneous emission may lead to a double 
peak structure in the spectrum; and (iii) we demonstrate 
how low-frequency noise leads to inhomogeneous broad- 
ening of the lasing peak. Finally, in Section IV, we study 
injection locking induced by an external coherent driv- 
ing, we discuss the main features of the spectrum, and 
we provide an estimate for the locking threshold. 



I. MODEL 

We consider a single-mode quantum resonator coupled 
to N a qubits (labelled by ji) and we account for both 
resonator and qubit dissipation. In the rotating wave 
approximation the system is described by the Hamilto- 
nian 

H = huj a^a+ -hujp^o-v + % VVcr^a + cr^a*) (1) 
+ (a + <J)X a + J2 {XW + X%o% + X^_) + H N . 

Apart from the photon operators, a and (V , we in- 
troduced the Pauli matrices acting on the single-qubit 



eigenstates a% 
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o~- = |0 M ) (1 M |. Dissipation is modelled by assuming that 
the oscillator and the qubits interact with noise opera- 
tors, X a and X£, X+, Xt, belonging to independent 
baths with Hamiltonian Hjsi in thermal equilibrium. The 
noise coupling longitudinally to the qubits, X^cr^, is re- 



sponsible for the qubits' pure dephasing. To describe 
lasing, we assume that a population inversion has been 
created in the qubits, which we describe by assuming that 
the effective temperature of the qubit baths is negative. 
In this way (for N a — 1) we model the essential prop- 
erties of the SSET laser used by Astafiev et al [8, 14]. 
Possible deviations from the standard Jaynes-Cummings 
oscillator-qubit coupling used in Eq. (1) are discussed 
in Appendix A. From the Hamiltonian (1), following the 
route described, e.g., in Ref. [16], we derive a set of quan- 
tum Langevin equations of motion, 



<n>/N 



2ig « a - aV) - ri(o£ - D ) + F?(t) , 
^+iuj )a-igY / V-+Fa(t). (2) 



dt 
d_ 

d 
dt.' 



Here the rate r x = r_i + I\ is the sum of excita- 
tion and relaxation rates, while Y v — T\/2 + T* is 
the total dephasing rate, which also accounts for pure 
dephasing due to the longitudinal noise described by 
r* In contrast to relaxation and excitation processes, 
Y* accounts for processes with no energy exchange be- 
tween qubit and environment, which thus do not af- 
fect the populations of the two qubit states. Further- 
more, k is the bare resonator damping. The param- 
eter Dq = (r-j- — Ti)/Ti denotes the stationary qubit 
magnetization in the absence of the resonator. In the 
present case, since we assume a negative temperature 
of the qubits baths and a population inversion we have 
Do > 0. The Langevin operators F?(t) with i = +, — ,z 
have vanishing averages and are characterized by their 
correlation functions, (Ff (t)Ff '(*')) = 5^D^g q (t - t'). 
The function g q (t — t') is assumed to decay on a time 
scale much shorter than the relaxation and decoherence 
times of the qubits and the oscillator. The diffusion co- 
efficients Df; are related to the rates introduced above, 

r T + r* (i + (<%)), Dt + = r i + r* (1 - (<#», 



D%_ 



d% z = 2T! - 2(r T - rj^}, d» + = 2r x «), d£_ = 

—2Y^(a l t). Similarly, the Langevin force F a (t) acting on 
the resonator with (F^(t)F a (t')) = KN^g a {t — t') is char- 
acterized by the rate K and thermal photon number -/V t h- 
The qubit and oscillator noises are assumed to be inde- 
pendent. For a further discussion of the different rates 
and diffusion coefficients we refer to standard quantum 
optics textbooks, e.g. Ref. [17]. 



II. STATIC PROPERTIES 

By taking appropriate products of Eqs. (2) and per- 
forming averaging, we arrive at the following equations 
for the average photon number, (n), the qubit polariza- 
tion (<j%) and the product (cr+a): 

±{a£) = -2z ff («a> - (a^ot)) - I^K) - D ) (3) 
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Figure 1: Color online - Scaled photon number, (n)/N a , in the 
threshold region versus scaled coupling, gs/N a for different 
values of N a . The other parameters are: u> p — u>o, Yi/ujo = 
0.016, r;/a>o = 0.004, Do = 0.975, k/lu = 3 ■ 10 -4 and 
Nth = 0. 
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u>o and the 
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where we introduced the detuning A 
total dephasing rate 7 = Y v + ~. 

In the stationary limit, neglecting the correlations be- 
tween different qubits, i.e. assuming (er^cr^) ~ (5 M „(1 + 
(0%)), the previous equations yield the following two ex- 
act relations between three quantities, the average polar- 



ization (S z (t)} with S z 
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If one of them is known, e.g., from a numerical solution 
of the Master equation, the other two can readily be de- 
termined. 

Factorizing the correlator, (S z n) m (S z )(n), on the 
right-hand side of Eqs. (6) gives results known in quan- 
tum optics as "semi-quantum model" [18]. It includes 
spontaneous emission processes, described by the term 
proportional to ((S z ) + 1). 

Spontaneous emission has a twofold importance for the 
issues described in the present work. First, at low tem- 
peratures it is responsible for the line- width of the lasers. 
Second, as noticed in Ref. [19], due to the low photon 
number, spontaneous emission is especially relevant in 
the description of the dynamics of single atom lasers. 
To illustrate this fact, we plot in Fig. 1 the scaled pho- 
ton number (n)/N a as a function of the scaled coupling 
g\/N a for different values of N a . The product g^/Na is 



kept constant to have a universal asymptotic behaviour. 
In the limit of large N a we observe a sharp lasing tran- 
si tion occur ring at the threshold coupling gthry/N a — 
y/n'y/(2Do), as predicted by the semiclassical theory. On 
the other hand, for low values of N a , and in particular 
for N a = 1, we find a smooth crossover between the nor- 
mal and the lasing regimes, which is due to spontaneous 
emission. In this case, although we cannot define a sharp 
threshold condition, we can still identify a transition re- 
gion centered at the semiclassical threshold coupling. 
The results presented in Fig. 1 were obtained analyti- 
cally using the semi-quantum approximation. In the case 
N a = 1 we compared such results with the numerical so- 
lution of the Master Equation and we obtained an agree- 
ment better than 10~ 3 . 



III. PHASE DIFFUSION 

For typical circuit QED parameters, i.e., for strong 
coupling g, the semi-quantum approximation, in spite of 
giving, as explained above, a good estimate of the sta- 
tionary photon number, cannot be used to study spectral 
functions. For the analysis of phase diffusion we thus pro- 
ceed with a hybrid approach: starting from Hciscnbcrg 
equations of motion we derive analytical expressions for 
the phase correlation time and frequency shift of the las- 
ing peak, both expressed as functions of the single-time 
averages, i.e., the photon number and qubit inversion in 
the stationary state. We then use the Master equation for 
the reduced qubit-resonator density matrix to calculate 
the single-time averages. 

From now on we consider a single-qubit laser, N a = 1, 
and analyze the laser and cross correlation functions 



0{t) = lim t ^ 00 (at(t + r)a(i)), 
G{t) = Umt^ooia+tt + rMt)). 



(7) 



Starting from the quantum Langevin equations (2) we 
derive a hierarchy of equations involving 0(t) and G(t). 
To truncate the hierarchy we split a(t) into an amplitude 
and phase, a{t) = \Jn(t) + 1 e~ tv ^\ and assume that 
the correlation time of phase fluctuations, 1/ Kd, is much 
longer than that of amplitude fluctuations, ~ 1/k [16]. 
This allows us to approximate for sufficiently long times, 

T > 1/K, 

Mt + r)ot(t + r)o(t)) * { -^^-(aHt + T)a(t)), (8) 

(Vn) 

while the correlator (a z y/n) can be estimated as 



{a z s/n) 1 / {<■>■■_>)) 



<V^> 2 



(9) 



Above threshold this approximation, which neglects 
terms of order (<5n 2 )/(n) 2 , is justified when the fluctu- 
ations of the photon number are much smaller than the 



average. Starting from Eqs. (2) and using the factoriza- 
tion (8) we obtain a coupled set of equations, 



J- t O{t) = (t«o-f)0(r)+* 5 G(r) 



<°w5), 



±G(t) = (i Wp -T v )G(T)-ig^^O(r). 



(10) 



These equations depend implicitly on all the parameters 
specified above. We focus on the case where the oscil- 
lator damping is much weaker than qubit's dephasing, 
k/2 <C r^, which is usually satisfied in singlc-qubit las- 
ing experiments. In this case we obtain from Eqs. (10) for 
the oscillator's spectral function, i.e., the real part of the 
Fourier transform of the correlator O(t), a Lorentzian, 
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and the frequency shift 

A 



SuJn = 
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The spectral function 0(u>) is proportional to the spec- 
trum measured in Ref. [8]. Equations (12) and (13) 
are the main results of the present work. Upon factor- 
ization of the correlator (cr z n) far above and below the 
lasing transition they reduce to results known from quan- 
tum optics [16]. The phase diffusion rate (12) is the sum 
of three terms. The first represents a thermal contribu- 
tion to the lincwidth and is negligible in the experimental 
regime explored in Ref. [8] . The second describes the ef- 
fect of relaxation processes of the qubit, which, due to the 
strong coupling, strongly increases the linewidth. The 
third term, due to quantum correlations, is essentially a 
measure of the coherent coupling between the qubit and 
the oscillator and leads to a reduction of the linewidth. 
Its effect is illustrated in the inset of Fig. 2, where we 
plot the diffusion constant Kd given by (12) covering the 
whole range from below to above threshold and compare 
it to the diffusion rate, k^ c , obtained by factorizing the 
correlator (a z n). As one can see correlations give signifi- 
cant quantitative corrections. Furthermore, we note that 
there is a wide range of parameters in which our approxi- 
mations remain consistent; in the case of strong coupling 
g this corresponds indeed to have Kd < K (see discussion 
around Eq. (8)). 

Fig. 2 also displays the dependence of the phase dif- 
fusion rate on the pure Markovian dephasing rate F* 
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Figure 2: Color online - Diffusion constant versus qubit- 
oscillator coupling and pure dephasing; parameters as in 
Fig.l. Inset: phase diffusion constants, Kd and K,f c , calcu- 
lated with and without taking into account correlation and 
the average photon number obtained from the Master equa- 
tion for F* v /lo = 0.004. 



showing a reduction of the linewidth with increasing T* 
above threshold. This surprising feature is a consequence 
of the fact that pure dephasing processes are not asso- 
ciated with emission of incoherent photons in the res- 
onator, and their main effect is simply a decrease of the 
effective qubit-resonator coupling. As one notes from the 
inset of Fig. 2, in single-qubit lasers far above threshold a 
reduction of the coupling has little effect on the saturated 
photon number but leads to a decrease of the incoherent 
photon emission rate, thus diminishing the linewidth. 

Another interesting manifestation of this effect is 
shown in the left panel of Fig. 3. Here we plot the height 
of the spectral line O(u>o + Sljq) = (n)/Kd as a function 
of the detuning, A. One can see that the optimal lasing 
conditions are realized somewhat out of resonance, 
where the effective coupling is weaker. We thus observe 
two peaks in the spectrum symmetrically shifted with 
respect to A = 0. Due to the strong coupling the photon 
number is roughly constant in the region between the 
two peaks. A similar structure in the output spectrum 
of single atom laser was also found in a numerical 
study by Ginzel et al. [20]. One might conjecture that 
this effect is the origin of the two spots observed in 
the experiment. It would explain why the peaks do 
not occur at resonance, but we have not succeeded to 
fit the experimental data in a satisfactory way [24]. 
In Ref. [8] it was proposed that the second peak is 
related to two-photon processes. Indeed deviations from 
the model used in Eq. (1), lead to an effective two 
photon-coupling between the qubit and the resonator. 
However, as described in the appendix, the two-photon 



coupling constant seems to be too small to produce any 
"two-photon lasing" . This results is also confirmed by 
the numerical solution of the Master Equation. We also 
note that in the experiment the second peak appears 
substantially shifted from the two-photon resonance 
condition, using the data of Ref. [8] , at the second peak 
we have 2u>q — uj p ~ OAujq. 

The linewidth of order of 0.3MHz observed in Ref. [8] 
is about one order of magnitude larger than what fol- 
lows from Eq. (12) (of the order of the Schawlow-Townes 
linewidth). Moreover in the experiment the laser line 
shows a Gaussian rather than a Lorentzian shape. Both 
discrepancies may be explained if we note that the qubit 's 
dephasing is mostly due to low- frequency charge noise, 
which cannot be treated within the Markov approxima- 
tion used in the derivation of Eqs.(ll)-(13). However, 
low-frequency (quasi-static) noise can be taken into ac- 
count by averaging the Lorentzian line in Eq. (11) over 
different detunings [22]. Assuming that the detuning 
fluctuations are Gaussian distributed with mean A and 
width a, such that r"i > a 3> Kd, we can neglect in the 
saturated limit the fluctuations of Kd and (n) and assume 
that the frequency shift Sujq depends linearly on the de- 
tuning A. From Eq. (13) we then have Scuq ~ Ak/(2T v ), 
and we obtain a Gaussian line with width a ~ <jk/(2T v ), 
where we remark that T v is the total markovian dephas- 
ing rate. The linewidth observed in the experiment is 
then reproduced by a reasonable value of a of order of 
300 MHz. In the case in which a is larger than Fi, 
the previous formula overestimate the linewidth since it 
doesn't take into account the decay of (n) out of thresh- 
old. In this case one can perform the averaging numer- 
ically. Anyway, in the presence of low-frequency noise, 
the linewidth is governed not by Kd as one may have ex- 
pected, but by Sujq. 



IV. INJECTION LOCKING 

We next investigate the behavior of the single-qubit 
laser when the oscillator is driven by an external laser 
field or seed light. To describe monochromatic driving 
with frequency Wdr and amplitude Eq 7 we add a term 
E ae luJdlt + H.c. to the Hamiltonian (1). It leaves the 
equations for the qubit operators unchanged, but modi- 
fies the quantum Langevin equation for the resonator, 



dt 



a = — I iloq 



a-t^e - *"*'*- 



■iga-+F a (t). (14) 



The average (a) now acquires a non-vanishing value and 
oscillates with the driving frequency u>di ■ In the resonant 



case Wdr 



too 



u) v , we estimate 



(a 1 ) 



(ae 



+iui dl t\ 



-iEo/Kd, where n d = \ (n d + \J{n 2 d + 2\E \ 2 /(n)) and 
Kd has the same functional dependence on (n) and (a z ) 
as in the undriven resonant case, Kd = § ( 1 + j3- 





Figure 3: Color online - Left panel: Maximum spectrum's 
amplitude O(u>o + Su>o) for different values of g. Right panel: 
Logarithm of the normalized resonator spectrum, \og[cjoO(uj)] 
in the presence of an external driving at resonance as a func- 
tion of the driving power, uq = 10GHz other parameters as 
in Fig. 1. 
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Finally, we consider the emission spectrum of the single 
qubit laser in the injection locking regime. For simplicity 
we neglect the low- frequency noise. In the double reso- 
nance regime, that is for Wdr = ^o = w p , we get a simple 
analytical expressions for 0(u>), 



0(w) = 



(U) - UJq) 2 + K 2 , 



■2%S(uJ-u>o)(a1)(a). (15) 



For low driving amplitude Eq, the resonator output is 
thus the superposition of two signals: the Lorentzian las- 
ing peak and a coherent peak due to the driving propor- 
tional to | (a) | . As one can see by combining the previous 
equation with the expression of (as), with increasing Eq, 
the height of the Lorentzian decreases and approaches 
zero, while the coherent peak grows. Eventually for large 
values of Eq, only the latter, which is amplified due to 
the coupling to the qubit, is visible in the spectrum. The 
driving amplitude, E , at which the Lorentzian peak dis- 
appears can be evaluated, using Eq.(15), as Eq ~ (n)R 2 d . 
The comparison to the experimental results of Ref. [8] 
can be done estimating the mean power absorbed by the 
system, Pg, as follows: Pg — E^ui^/Rd [17]. To illustrate 
the locking transition we plot in Fig. 3 (right panel) 
the spectrum (15). The output spectrum is centered at 
u> = since 8ujq — 0. In the numerics we assumed that 
the injected signal has a Lorentzian shape with width 
Wdr = ft/200. In the presence of a finite qubit-oscillator 
detuning A = uj p — loq, the position of the lasing peak 
would be shifted. 



V. CONCLUSION 

We analyzed in detail the spectral properties of single- 
qubit lasers. Our main conclusions are: 
- Due to the strong coupling nontrivial structures appear 
in the spectrum, which are not visible in the average pho- 
ton number. As shown in Fig. 3 the optimal lasing con- 
ditions are realized for two values of the detuning, which 



are symmetrically shifted from A = 0. At these two hot- 
spots the output spectrum has the maximum height and 
it is centered around the frequency uiq ± 8u)q. 

- Low-frequency noise strongly affects the line shape of 
the two peaks, leading to an inhomogeneous broadening. 
In comparison, the natural laser lincwidth due to spon- 
taneous emission is negligible. 

- Although we did not produce a quantitative fit to the 
data, we presented a possible explanation of the double- 
peak structure observed in the experiment of Ref. [8]. 
We obtained an estimate of the lincwidth due to low fre- 
quency noise in qualitative agreement with the exper- 
iments, and we evaluated the locking threshold in the 
injection locking experiment. 

We acknowledge fruitful discussions with O. Astafiev, 
A. Fedorov and M. Marthaler. The work is part of the 
EU 1ST Project EUROSQIP 



Appendix A 

Here we briefly discuss the validity of the model intro- 
duced in Section I, when applied to describe the experi- 
ment of Astafiev et al. [8] . The single-qubit laser realized 
in Ref. [8] , consists of a charge qubit coupled capacitively 
to a single-mode electrical resonator and can be thus de- 
scribed in the qubit 's eigenbasis by the Hamiltonian 



1 



,n 



H = —AEa z + hu>oa'a—hgo(sm(a z +cos(a x ) [a -+ 

(Al) 
The angle £ and the qubit energy splitting depend on the 
charging and Josephson energies, e c h and Ej, tan £ = Zjjp- 

and AE = \/e 2 ch + E 2 . In order to identify the one- and 
two-photon couplings, we now apply a Schrieffer- Wolff 
transformation U = e lS with S = i 2si ^<Jz (a - at) 
and perform a perturbation expansion in the parameter 
go/ojQ. The transformed Hamiltonian, H = U' HU, thus 
becomes 

H ~ -AEo z + huj a^a + hg x a x (a + a f ) (A2) 



+hg2icr y (a 2 — (a^)' 



Here we neglected terms of order (go/u)o) and intro- 
duced the two coupling constants g\ = —go cos £ and 

2 2 

.92 = — -^ sin C cos £ for one-photon and two-photon 
transitions, respectively. For the parameters used in the 
experiment the coupling gi is roughly two orders of mag- 
nitude smaller than the one photon coupling and be- 
low the s emiclassical th reshold for the two-photon lasing, 
92 = V^^tp/fiiDo) [23] . In the parameters regime ex- 
plored in the experiments, the Hamiltonian used in Eq. 
(1) gives thus a good description of the dynamics of the 
system. 
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